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We predict the existence of spatiotemporal solitons ( "light 
bullets" ) in two-dimensional self- induced transparency media 
embedded in a Bragg grating. The "bullets" are found in an 
approximate analytical form, their stability being confirmed 
by direct simulations. These findings suggest new possibilities 
for signal transmission control and self-trapping of light. 
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Light propagation in periodic dielectric structures ex- 
hibits a variety of unique regimes that are technologicaUy 
promising: nonlinear fihering, switching and distributed 
feedback amplification ||l[ . Of particular interest are gap 
solitons (GSs), i.e., moving or standing self-localized field 
structures centered in a band gap of the grating. These 
self-localized field structures arise due to the interplay 
between the medium nonlinearity and resonant Bragg re- 
flections. Their spectrum is tuned away from the Bragg 
resonance by the nonlinearity at sufficiently high field in- 
tensities. Theoretical studies of gap solitons in Bragg 
gratings with Kerr nonlinearity |^ have been followed up 
by their experimental observation in a nonlinear optical 
fiber with the grating written on it P] . Gap solitons have 
also been theoretically studied in gratings with second 
harmonic generation B. 

A principally different mechanism giving rise to gap 
solitons has been found in studies of models consisting of 
a periodic array of thin layers of resonant two-level sys- 
tems (TLS's) separated by half-wavelength nonabsorbing 
dielectric layers, i.e., a resonantly absorbing Bragg reflec- 
tor (RABR) §-§ (Fig. |I|). Such a RABR has been 
shown, for any Bragg reflectivity, to produce a vast fam- 
ily of stable solitons, both standing and moving ones. As 
opposed to the 27r solitons arising in self-induced trans- 
parency (SIT), i.e., near-resonant field-TLS interaction 
in a uniform medium [g-10|, gap solitons in a RABR 



can have an arbitrary pulse area |g|-|7|]. The existence of 
GS solutions can only be consistently demonstrated in a 
RABR with thin active TLS layers. By contrast, a re- 
cent attempt |1^ to obtain such solutions in a periodic 
structure uniformly filled with active TLS's is physically 
doubtful, and fails for many parameter values. 

The potential applications of GSs are based on the sys- 
tem's ability to filter out (by Bragg refiection) all pulses 
except for those satisfying the GS dispersion condition, 



as well as the control of pulse shape and velocity [^-0. 
It would be clearly desirable to supplement these advan- 
tages by immunity to transverse diffraction of the pulse, 
i.e. to achieve simultaneous transverse as well as longitu- 
dinal self-localization in the structure. This calls for the 
consideration of "light bullets" (LBs), multi-dimensional 
solitons that are localized in both space and time. In the 
last decade they have been theoretically investigated in 
various nonlinear optical media |l2|-fi8| , and the first ex- 
perimental observation of a quasi two-dimensional (2D) 
bullet was recently reported [^. In a recent work P0[ , 
we have predicted that uniformly 2D and 3D self-induced 
transparency (SIT) media can support stable light bul- 
lets. 

In the present work, we aim to extend this investigation 
to RABR's of the kind shown in Fig. || so as to combine 
LB and GS properties in resonantly absorbing media. 
We find that a RABR with any Bragg refiectivity and 
any absorption cross-section can support the propagation 
of attenuated stable LBs, which are closely related to 
the bullets we have found in uniform SIT media EG]. It 
should be noted that 2D LBs supported by a combination 
of the Bragg refiector with a different (second-harmonic- 
generation) nonlinearity were theoretically investigated 
in Ref. |l|]. 

We start by considering a 2D SIT medium with a 
spatially-varying refractive index n{z,x), which is de- 
scribed by the lossless Maxwell-Bloch equations M , 



i£^^ + n^£r+£z+i{l-n^)S -V = 0, 

Vr-£W ^ 0, 
1 



Wr 



-{E*V + V*E) = 0. 



(la) 
(lb) 

(Ic) 



Here £ and V are the slowly varying amplitudes of the 
electric field and polarization, W is the inversion, z and 
X are the longitudinal and transverse coordinates (mea- 
sured in units of the resonant-absorption length), and 
T is time (measured in units of the input pulse dura- 
tion). The Fresnel number, which governs the trans- 
verse diffraction in the 2D and 3D propagation, has 
been incorporated in x and the detuning of the car- 
rier frequency wq from the central atomic-resonance fre- 
quency was absorbed in £ and V . The Fresnel num- 
ber F, detuning Afi, and wave vector /cq = cjo/c can 
be brought back into Eqs. (1) by the transformation 
£{t, z, x) — > (2/fco) £{t, z, x) exp(— iAil • r), Vlr, z, x) -^ 



(4//c2) V{t, z, x) exp( 



lArJr), W{t,z,x) 
(2/fco) z, and x 



{'i/kl) W, 



T -^ (2/fco) T,z^ (2/fco) z, and x -^ ^2/koFx. To ne- 
glect the polarization dephasing and inversion decay we 
assume pulse durations that are short on the time scale 
of the relaxation processes. Equations (1) are then com- 
patible with the local constraint [T'P -I- W^ = 1, which 
corresponds to the Bloch- vector conservation |^. In the 
absence of the x-dependence and for n{z, x) = 1, Eq. ( p^ ) 
reduces to the sine-Gordon (SG) equation which has the 
soliton solution £(t, z) = 2a sech(ar — z/a + Go), where 
a and Qq are real parameters. 
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FIG. 1. Periodic modulation (g) of the refractive index in 
a structure of periodically alternating layers. Adapted from 
Ref §. 

We proceed to search for LBs in a 2D medium subject 
to a resonant periodic longitudinal modulation, i.e., the 
above-mentioned RABR (Fig. p. Accordingly, we as- 
sume the following periodic modulation of the refractive 
index [ill 



^(z) — n^ [1 + ai cos(2fcc2;)] . 



Here, no and ai are constants and fcc = iOc/c, with ujc 
being the central frequency of the band gap. A RABR 
is then constructed by placing very thin layers (much 
thinner than 1/fcc) of two-level atoms, whose resonance 
frequency is close to Wc, at the maxima of this modu- 
lated refractive index. We aim to consider the propaga- 
tion of an electromagnetic wave with a frequency close 
to ijOc through a 2D RABR. Due to the Bragg reflec- 
tions, the electric field S is decomposed into forward- and 
backward-propagating components £p and £s, which 
satisfy equations that are a straightforward generaliza- 
tion of the ID version in Refs. 
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-F7?S+^ + TfT.+ - 2Vr - 2ir]V 


= 0, 


(3a) 


~^^rxx + ^'^zxx + ^TT ^ ^zz 






~ri^-, + 7fi:- +2Vz 


= 0, 


(3b) 


Vr + ilS.VLV - S+VF 


= 0, 


(3c) 


Vi/^ + i(S+*-p + S+7?*) 


= 0. 


(3d) 



Here Yr^ = 2To/ino(f_F ± ^sM'h, tq = uq^l '^yjhplmj^, 
H is the transition dipole moment, and po is the density 



of the two-level atoms. The parameter 77 is a ratio of 
resonant-absorption length in the two-level medium to 
the Bragg reflection length, and can be expressed as i] — 
aiWcTo/4. 

In the ID case, a family of exact soliton solutions to 



was 



Eqs _ 

found in Ref. @: S(±) == 2 (a; -Va^ -I- 2) seche(r, z) 
exp(i$), V = 2scch6(r, 2) tanh9(T, z) • exp_(j$), W = 
2sech^9(T, z) — 1, where 6(t, z) = ar + \/ o? + 2z + 9o, 
$ = rjMr + rjNz + (J), M = -(a^ + 1), N = -a^Ja^ + 2, 
and Afi = ^{a^ + 1). The shape of the fields 5]+ and Yr 
in this solution is similar to the SG soliton in the uni- 
form ID SIT medium. Inspired by this analogy and the 
fact that there exist LBs in the uniform 2D SIT medium 
which reduces to the SG soliton in ID ||22l, we search for 
a LB solution to the 2D equations (g), which reduces to 
the exact soliton in ID. To this end, we try the following 
approximat ion , 

S+ = 2aVsecheiseche2 e^■^^Mr+^^Nz+^^|^ ^ 

Yr 



-2 \Jo?- + 2 VsechGiseche 



2 e 



ir) M T -{-ir] N z-\'i7r / 4 



V = Vscche 1 scche2 { (tanh 6 1 +tanhe2)^ + 
-a^C"' [(tanh 61 - tanh 62)^ - 2(sech^ei 



sech^Ga)]'} 



1/2 iriMT+ir)Nz+iu 



w = [i-|-pp] 



1/2 



with 01 (t, z) 



ar + vc? 
f Oo-Ce, 



+ 2z + eo + Cx, e2(T,z) = 

the phase v and coefficient C 



ar + \J o? + 2 z 
being real constants. 

The ansatzM) satisfies Eqs. (pa) and (3b) exactly, 
while Eqs. (|3d|) are satisfied to order |a|C, which re- 
(2) quires that \a\C^ ^ 1. The ansatz is relevant for arbi- 



trary 77, admitting both weak (77 <C 1) and strong (77 > 1) 
reflectivities of the Bragg grating, provided that the de- 
tuning remains small with respect to the gap frequency, 
or ?7 ^ bJcjio? + 1). Gomparison with numerical simu- 
lations of Eqs. (^, using (^ as an initial configuration, 
tests this analytical approximation and shows that it is 
indeed fairly close to a numerically exact solution; in par- 
ticular, the shape of the bullet remains within 98% of 
its originally presumed shape after having propagated a 
large distance, as is shown in Fig. ||. 
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FIG. 2. The forward-propagating electric field of the 2D 
"hght buUet" in the Bragg reflector, \£f\, vs. time r and 
transverse coordinate x, after having propagated the dis- 
tance z = 1000. The parameters are a = 1, C — 0.1 and 
00 = —1000. The field is scaled by the constant h/4:TofJ,no. 

We stress that 2D or 3D LB solutions of the variable- 
separated form S+ ^ S^ '^ /(r, z) ■ g{x) do not ex- 
ist in RABR. Indeed, substituting this into Eqs. (Ba) 



and (3b) yields only a trivial solution of the form E ~ 
exp {iAt + iBx), with constant A and B. 

We briefly discuss experimental conditions under 
which LBs can be observed in RABRs. The incident 
pulse should have uniform transverse intensity within its 
diameter d. For the transverse diffraction to be strong 
enough, one needs aeScP/^a < 1, where aos and Xq 
are the inverse resonant-absorption length and carrier 
wavelength, respectively p^. For acff ~ lO'^ m~^ and 
Ao ^ 10~^ m, one thus requires a diameter d < 10~^ 
m, which implies that the transverse medium size La; ~ a 
few /im. 

In order to realize a RABR, thin layers of rare-earth 
ions |23| embedded in a semiconductor structure with 
a spatially-periodic RI |^J] may be used. The two-level 
atoms in the layers should be rare-earth-ions with density 
of 10^^ — 10^^ cm^'^, whose resonant-absorption time and 
inverse length are respectively tq ~ 10""'^^ — 10~^^ s and 
ftoff ~ 10* — 10^ m^^. The parameter 77 can vary from 
to 100 and the detuning is ^ 10^^ - 10^^ s"^ In a 
RABR with transverse size 10 /im, LBs depicted in Fig. g 
are localized on the time and transverse-length scales ^ 
10^^^ s and 1 /im. Cryogenic conditions strongly extend 
the dephasing time T2 and thus the LB lifetime, well 
into the /isec range [ESl. The construction of suitable 
structures constitutes a feasible experimental challenge. 

To conclude, we have studied light bullets in SIT me- 
dia embedded in a resonantly-absorbing Bragg reflector. 
Light bullets in a multi-dimensional Bragg reflector have 
the potential of serving as a novel type of optical filters, 
which stably transmit selected signal frequencies through 
their spectral gap and block others. They can also be 
used to both spatially and temporarily localize light in 
certain frequency bands. 
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